Exact solutions of an »S0(5)-invariant spin-3/2 Fermi gas model 
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An exactly solvable model describing the dilute spin-3/2 fermion gas in one-dimensional optical 
trap is proposed. The diagonalization of the model Hamiltonian is derived by means of the Bet he 
ansatz method. Exotic spin excitations such as the heavy spinon with fractional spin 3/2, the neutral 
spinon with spin zero and the dressed spinon with spin 1/2 are found based on the exact solution. 
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The study on the optically and magnetically trapped 
ultracold atoms has attracted a lot of attentions in the 
recent years. By using the magnetic fields or laser beams, 
the atoms can be trapped and cooled down to very low 
temperatures 0, 0, 0, 0, 0, 0, 0, 0, 0, Ej • A fascinating 
fact is that not only the components but also the physical 
parameters of such systems can be manipulated in exper- 
iments, which allow us to mimic the ordinary correlated 
electron systems and even to explore new states of mat- 
ter. For examples, by means of the Feshbach resonance 
techniques, one can tune the scattering lengths or the 
interactions among the atoms; putting an optical lattice 
on one can study the behaviors of the atoms in a tunable 
periodic potential and with strong anisotropic traps, the 
one- and two-dimensional quantum systems can be real- 
ized. One of the hot research topics in this field is the 
study on the high spin cold atom systems. The physics 
of bosonic ultracold atoms with hyper fine spi n F = 1 
have been extensively studied [HI, [lH, EH, (3, EE 
The degenerate Fermi gas with hyperfine spin F = 3/2 
could be obtained by cooling alkali atoms 132 Cs, as well 
as alkaline-earth atoms 9 Be, 135 Ba and 137 Ba in exper- 
iments. The ground states of these high spin fermionic 
atoms in optical traps are very rich. Many peculiar quan- 
tum orders and exotic collective excitations appear in 
these systems, which are rare in the ordinary interacting 
electron systems [I?}- Several methods have been devel- 
oped to approach these ultracold atom systems. For ex- 
amples, Wu et al. showed that the one-dimensional (ID) 
spin-3/2 Fermi gas with s-wave scattering possesses the 
50(5)-invariance and obtained the phase diagram of the 
system by means of bosonization [18|, 19] : the functional 
integral approach was applied on the interacting spin-3/2 
fermionic ultracold atoms in 2D square optical lattice to 
study a variety of Mott insulating phases 0, EH; two dif- 
ferent superfluid phases in F = N — 1/2 fermionic cold 
atom system were found, i.e., an unconfined BCS pairing 
phase and a confined molecular-superfluid made of 2N 
fermions, depending on whether a discrete symmetry is 
spontaneously broken or not [22]. 

Despite the fast progress in this field, the understand- 
ing to the high-spin correlated systems is still far from 
satisfied. Reference models with exact solutions are un- 



doubtedly needed to get deep insight about these in- 
teresting quantum many-body systems. An important 
progress in this aspect is made by Controzzil and Tsve- 
lik, who constructed an exactly solvable model for isospin 
F — 3/2 fermionic system with linear dispersion relations 
[23I ]. In this Letter, we propose a new exactly solvable 
model which describes properly a ID spin-3/2 cold atom 
system. With the Bethe ansatz solutions, we derive out 
exactly the ground state. Exotic spin excitations such as 
the heavy spinons with fractional spin 3/2, the neutral 
spinons with zero spin and the dressed spinons with spin 
1/2 are found. 

The fermions with hyperfine spin F = 3/2 in a ID 
trap is appropriately described by the following model 
Hamiltonian 



N 



3=1 



d 2 



N 



gmSixj-x^, (1) 



l<j,m 



where V e (xj) is the external potential, g m is the two- 
body coupling constant in the total spin m = 0, 2 chan- 
nel and is the projection operator onto the spin m 
channel. Because of the antisymmetric wave functions of 
the fermions, non-trivial scattering processes occur only 
in the m = and m — 2 channels, and those in the m = 1 
and m = 3 channels are irrelevant for the (5-function in- 
teraction potential. In the low density case, we omit the 
external field. The Hamiltonian (pQ) can be rewritten as 
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where cq = go/3 + 2g 2 /3, c 2 = g 2 /3 — go/ 3 and Sj is the 
spin-3/2 operator. If c 2 = 0, the system (j2j) is SU(4:)- 
invariant and was solved exactly by Sutherland [24(. In 
the present Letter, we show that the system ([2]) pos- 
sesses another integrable line cq = c/2, c 2 = — 2c/3 at 
which the physical properties are quite different to those 
of the Sutherland model [24], [25[ . An obvious fact is that 
in our case, the particle number of an individual spin 
component is no longer conserved because of the broken 
6'/7(4)-symmetry. Nevertheless, the system (j2j) is still 
6'0(5)-invariant at this new integrable line, and the fol- 
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lowing three independent conserved quantities are hold: following form in the thermodynamic limit 
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where N s indicates the particle number with the spin 
component s = ±1/2, ±3/2. 

Assume the wave function takes the following form 

*(xisi,--- ,x N s N ) = ^20(x Ql < ••• < x Qn ) 

xAs^sMPyZiLrk^xo^ (4) 

where Q = (Qi, • • • , Qn) and P = (Pi, • • • , P/v) are the 
permutations of the integers 1, • • • , N; kj are the quasi 
momenta carried by the particles; 0(xq 1 < • • • < xq n ) = 
Q( X Q N - x Qn-i) ' ' ' Q( X Q2 - x Qi ) and 0(x - y) is the step 
function. With the standard coordinate Bethe ansatz 
method, we obtain the two-body scattering matrix as 
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which satisfies the Yang-Baxter equation 

Si2(*i - k 2 )S 13 (h - k 3 )S 23 (k 2 - k 3 ) 

= S 23 (k 2 - k 3 )S 13 (h - k 3 )S 12 (h - k 2 ). (6) 
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where A^ and pi^ are the real parts of the n-string 
of A and the m-string of /i, respectively. Denote p, 
r] n and a m as the densities of A n-strings and p m- 
strings in the thermodynamic limit N — > oo, L —> oo 
and N/L — > finite, respectively, and p h , r\\ and the 
corresponding densities of holes. At temperature T, the 
Gibbs free energy of the system with an external mag- 
netic field ft and chemical potential A is 



F = E — AN — ft [ -N - Mi — M 2 ) - TS, 



(10) 



where E = L J k 2 p(k)dk is the energy; N = L J p(k)dk 
is the particle number; Mi = L^2 n n J r] n (X)dX; M 2 = 

L Em m I VrnfaWl S = L f{(p + p h ) Hp + P h )-plnp- 
P h In p h + Y. n l(Vn + Vn) HVn + Vn) ~ Vn In Vn ~ Vn ^ Vn\ ± 

E m [K ± O ln(cr m + <j£j - cr m In a m - In a^]}dk 
denotes the entropy. Minimizing the Gibbs free energy 
at the thermal equilibrium, we obtain the following ther- 
modynamic Bethe ansatz equations (TBAE) 



Applying further the nested algebraic Bethe ansatz with 
periodic boundary conditions [26], we obtain the follow- 
ing Bethe ansatz equations(BAE): 
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where = N 3/2 + iV 1/2 ± A^_i /2 + A± 3/2 , Af x = iV 1/2 ± 
2A^_ 3/2 + N_ 1/2 , M 2 = N_ 1/2 + iV_ 3/2 and L is the 
length of the system. The corresponding eigen energy of 
the Hamiltonian is E — J^jLi kj- 

We consider the c > 0, i.e., the repulsive interaction 
case. By carefully checking the structure of the BAE 
we find that all the charge rapidities kj take real values, 
indicating the absence of charge bound state. However, 
the spin rapidities \p and p v may form strings with the 
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where p = p h /p, fj n = Vn/Vn, Vm = iK, G n 
a n /(ao + a 2n ), a = (5(x), a n (x) = 4nc/[7r(16x 2 + n 2 c 2 )] 
and a n * /(a;) = / a n (x - y)f(y)dy. 

The ground state configuration of the system can be 
obtained by taking the limit of T — » and ft — > 0. In this 
case, most of the string densities are zero and the TBAE 
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(fTTj) are reduced to 

P( k ) = ^~ + a l * Vl( k ) + a 2 * 7/2 (fc), 
Z7T 

771 (A) = / ai (A — k)p{k)dk — a 2 * 771(A) 

-(ai + a 3 ) * 772(A) + ai * cr(A), (12) 
rQ 

7/2(0 = / a 2(^ - k)p(k)dk - (ai + a 3 ) * 771(0 
«/-Q 

-(2a 2 + a 4 ) * 772(0 + (ai + ^3) * 0"(O> 
cr(/i) = ai * 771 (/x) + (ai + a 3 ) * 772 - a 4 * cr(/i). 

The Fermi point Q is determined by the density of par- 
ticles N/L = J Q Q p(k)dk. M 1 /L = /[771(A) + 2r7 2 (A)]dA 
and M 2 /L = J a(p)dp. Such a ground state configura- 
tion is quite different from that of the SU (4) Sutherland 
model, where there is no string or spin bound state in 
the ground state. In the present £0(5) case, part of 
the spin rapidities form 2-strings which heavily affect the 
spin excitations as we shall show below. From the solu- 
tions of Eq. ([T2]) . we can easily derive out Mi = N and 
M2 = AT/2, which give the total spin of the ground state 
S = 3N/2 — M1—M2 = 0, indicating a spin singlet ground 
state. When c tends to zero, the function a n (x) — > S(x). 
In this case we recover the free Fermi gas solutions 

= I*?, "ig. (is) 

7T L 07T L 7T 

When c — > 00, a n (x) — > 0, we recover the Tonks- 
Girardeau solutions 

Based on the ground state configuration, the elemen- 
tary excitations of the system can be studied exactly. 
In the integrable models, the excitation energies are 
uniquely defined by the so-called dressed energies e = 
Tin/?, Cn = Tln77 n and £ m = Tln<j m . In our case, when 
T — > 0, only e, ("1,2 and £1 are left. The dressed energies 
of the ground state for Q = 1 and c = 1 are shown in 
FigJTJ The low-lying excitations of the system can be 
studied systemically by adding particles, holes or strings 
into the ground state configuration of the rapidities. The 
excitation energy reads 

n h n p mi 

AE = £ e (^) + £ e (fcJ) + ^Ci(At-) 
j=l j=l j=l 

■mi m 3 

where n^, n p and mi 5 2,3 are the numbers of the charge 
holes, excited charges, holes in real A sea, in 2-string A 
sea and in real p sea, respectively; kj, fcj, A^ 2 j - and p!- 
are the positions of the corresponding charges and holes. 
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FIG. 1: (Color online) The dressed energies of the ground 
state (Q = 1, c = 1). e, (1, ("2 and £ are the dressed energies 
of rapidities fc, real A, 2-string A and real /i, respectively. At 
the Fermi point Q, the dressed energy e(Q) is zero. 




FIG. 2: (Color online) The low- lying excitations (c = 1, Q — 
1). AE and P are the energy and the momentum carried 
by the excitation, (a) The charge-hole excitations, (b) The 
excitation of two holes of real A. (c) The excitation of two 
holes of real p. (d) the excitation of four holes of A 2-string. 

Formally, the extra strings contribute nothing to the en- 
ergy because the contribution of such strings is exactly 
canceled by the rearrangement of the Fermi sea. Some 
of the low-lying excitations are shown in FigJSJ The real 
k charge-hole excitation in our case is similar to that of 
6'/7(4)-invariant Sutherland model as shown in FigJSfa). 
However, the excitations in the spin sector are quite dif- 
ferent from those of the Sutherland model. The spin 
quanta carried by the spin excitations reads 

S = -m 2 + 2m 3 + ^(2 - Z)ra A(0 + ^(1 - t)m^ t) , (16) 

l>3 t>2 

where m X (i) and m^t) are the numbers of A /-strings and 
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FIG. 3: (Color online) The single- hole excitations (Q = 1, 
c = 1). AE and P are the energy and the momentum carried 
by a single hole, respectively, and n is the density of particles. 
The solid line is the single- hole excitation of real k. The 
dashed line is that of real A. The dotted line is that of A 
2-string and the dotted dashed line is that of real ji. 



TABLE I: Possible hole configurations. 
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neutral spinons). The A 2-string hole pair can not ex- 
ist independently but must be associated with a neutral 
spinon (the second column in Table [Tj) or a real \i hole (the 
third column in Table [Tj). If we add further a A 4-string 
into a 2-string hole pair and a real A hole configuration, 
the total spin of this excitation is 1. In this case each of 
the A 2-string holes carries a spin 1/2. Such a dressed 
hole is quite similar to the ordinary spinon (named as 
dressed spinon). Four A 2-string holes may exist inde- 
pendently (the fourth column in Table [Tj) . If we put fur- 
ther one A 6-string and one \i 3-strings into this four hole 
configuration, we get the SO (5) spin singlet excitation. 
The simplest excitation in the \i sector is a pair of real 
fi holes (the fifth column in Table [Tj). This excitation is 
quite similar to a real A hole pair but each of the real fi 
hole carries a spin 2. Joint pair of a real A hole and a real 
H hole may also happen as shown in the last column in 
Table H Other kinds of spin excitations can be obtained 
by analyzing the BAE but most of the complex excita- 
tions are generally composed of the excitations listed in 
Table H 
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the fi t-strings formed in the excitations, respectively. 
We note the numbers of holes and strings added are not 
independent but satisfy some constraints determined by 
the BAE: 
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where AM} 2 (integers) indicate the number changes of 
A, \i /-strings. Several possible hole configurations are 
listed in Table H Nevertheless, the energies of the holes 
are additive as shown in Eq. ([T5j) and the thermodynamic 
behavior of the system is mainly determined by the dis- 
persion relations of the individual holes, which are shown 
in FiglSl Interestingly, the single A 2-string hole carries 
the lowest energy with spin 3/2 (named as heavy spinon) 
and therefore dominate the low temperature thermody- 
namics of the system. 

The simplest spin excitation is a real A hole-pair (the 
first column in Table [Tj) , corresponding to the two do- 
main walls of a single excited domain. However, unlike 
the usual spinons, such holes carry zero spin (named as 
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